Half-metallic ferromagnetism and spin-polarization in Cr02 
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We present electronic structure calculations in combination with local and non-local many-body 
correlation effects for the half- metallic ferromagnet Cr02. Both zero temperature results from a 
Variational Cluster Approach as well as finite-temperature Dynamical Mean Field Theory results 
support the existence of low-energy non-quasiparticle states, which turn out to be essential for a 
correct quantitative description of depolarisation effects. 



In recent years new magnetoelectronic devices have 
emerged that are able to exploit both the charge and 
the spin of electrons by means of spin-polarized currents 
and spin-dependent conduction [ij, t2] ■ The capabilities of 
these devices depend crucially upon the spin polarization 
of the underlying material. In this area half-metallic fer- 
romagnets are of particular interest, especially chromium 
dioxide Cr02 is a strong candidate for low-field magne- 
toresistance devices with a favorable switching behaviour 
at small fields 

Based on first-principle calculations Cr02 is classified 
as a half-metallic ferromagnet 4] . Point contact measure- 
ments at superconductor metal interfaces reveal a spin 
jolarization of the conduction electrons larger than 90% 
5i, ^S], supporting the half-metallic nature predicted by 
band theory Q. In addition, Cr02 fulfils yet another 
essential requirement for practical purposes, namely a 
high Curie temperature,^ determined experimentally in 
the range of 385-400K 8]. The saturation magnetiza- 
tion at lOK was reported to be 1.92/iB per Cr site 
which is close to the ideal value {2^.b) and is consistent 
with half-metallicity. The magnetic susceptibility in the 
paramagnetic phase shows a Curie- Weiss like behavior 
indicating the presence of local moments suggest- 
ing a mechanism of ferromagnetism beyond the standard 
band or Stoner-like model. 

The band structure of Cr02 has been calculated self- 
consistently within the local-spin-density-approximation 
(LSDA) by several authors. K.-H. Schwarz [7| first pre- 
dicted the half-metallic behavior for this material with 
a spin moment of 2/iB per formula unit. Later on, sev- 
eral experiments including photoemission 11 ill, soft x-ray 
absorption [l^, resistivity and optics [1^ have sug- 
gested that electron correlations are essential to under- 
stand the underlying physical picture in Cr02. In this 
spirit, LSDA-I-U calculations 'l5'| have been used to study 
the effect of large on-site Coulomb interactions. These 
studies suggested that Cr02 is a material with a negative 
charge transfer gap which leads to self-doping. In con- 
trast to the on-site strong correlation description, trans- 
port and optical properties obtained within the density 
functional theory using LSDA and generalized-gradient 



approximation [la|, suggest that the electron-magnon 
scattering is responsible for the renormalization of the 
one-electron bands. More recent model calculations also 



suggest the importance of orbital correlations jlTI. Il8|. 

An appropriate treatment of correlation effects, which 
goes beyond the mean-field LSDA-I-U treatment is 
achieved in the framework of Dynamical Mean Field The- 
ory (DMFT) J^] , in which a local energy-dependent self- 
energy is adopted. However, nonlocal self-energy effects 
turn out to be important for a correct description of the 
so-called non-quasiparticle (NQP) states ^20,,^], which 
are crucial for the depolarisation in half-metallic ferro- 
magnets. 22, |2^ [3 For this reason, we present 
here for the first time a realistic electronic-structure cal- 
culation for Cr02 , where non-local self-energy effects are 
treated within the Variational Cluster Approach (VCA) 
[H, S H lii- While the VCA calculation is cur- 
rently limited to zero temperature, we also present finite- 
temperature results obtained within a DMFT descrip- 
tion. 

In order to achieve a realistic description of the mate- 
rial, the parameters of the correlated model Hamiltonian 
are determined via ab-initio electronic-structure calcula- 
tions. Cr02 has a rutile structure with Cr ions forming 
a tetragonal body-center lattice. Cr^+ has a closed shell 
Ar core and two additional 3d electrons. The Cr ions are 
in the center of the CrOe octahedra. Therefore, the id 
orbitals are split into a t2g triplet and an excited eg dou- 
blet. With only two 3d electrons, important correlation 
effects take place essentially in the t2g orbitals, while eg 
states can be safely neglected. The cubic symmetry is 
further reduced to tetragonal due to a distortion of the 
octahedra, which partially lifts the degeneracy of the t2g 
orbitals into a dxy ground state and d^z+zy and dxz-zy 
excited states [15l |30|. We therefore restrict the VCA 
calculation to three orbitals on each site representing the 
Cv-t2g manifold described by the Hamiltonian 



H — Hq + — Umm'm"m"'cl„ 



i{m.(7} 
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where Cimcr destroys an electron with spin a on orbital 
m on site i. Here, Ho is the noninteracting part of the 
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Hamiltonian restricted to the t2g orbitals obtained by 
the downfolding procedure implemented within the N- 
order Muffin-Tin Orbital (NMTO) method using 
the Local Density Approximation (LDA). In this way, 
the effects of the remaining orbitals are included effec- 
tively by renormalisation of hopping and on-site energy 
parameters. The on-site Coulomb-interaction terms in 
Eq. Ill) are expressed in terms of the diagonal direct 
coupling Ummmm = U , the off-diagonal direct coupling 
Ujnjn'Tnm' = U' and the exchange coupling Umm'm'm = J 
[331 ■ Notice that spin-rotation invariance is automati- 
cally guaranteed by the form Eq. [H i. e. spin-flip terms 
are also included in our calculation. The pair-hopping 
term is also included via Ummm'm' — J ■ The Coulomb- 
interaction (C7) and Hund's exchange (J) parameters be- 
tween t2g electrons are evaluated from first principles by 
means of a constrained LSDA method [34]. While this 
method is not very accurate, and may produce slightly 
different values depending on its exact implementation, 
it is quite clear that for Cr02 the interaction parameters 
are not too strong. An appropriate choice is C/ « 3 eV 
and J « 0.9 eV. Corrections for a "double-counting" of 
the interaction only produces an irrelevant constant shift 
of the chemical potential [sll, as we are considering a 
model Hamiltonian with fixed number of electrons. 

The VCA accesses the physics of a lattice model in 
the thermodynamic limit by optimizing trial self-energies 
generated by a reference system. The reference system 
consists of an isolated cluster having the same interaction 
as the original lattice, but differing in the single-particle 
part of the Hamiltonian. An approximation to the self- 
energy of the lattice model is then obtained by finding the 
saddle point of an appropriate grand-canonical potential 

Q. = VL' + Tr In GvcA - Tr In Gd (2) 

with respect to the single-particle parameters of the ref- 
erence system. Here, Gd is the Green's function ma- 
trix of the reference system (cluster) , which is calculated 
numerically as a function of frequency and spin, using a 
zero-temperature Lanczos procedure. The lattice Green's 
function GycA is obtained by a matrix form of the Dyson 
equation, whereby the self-energy E of the reference sys- 
tem is used as an approximation to the lattice self-energy. 
Finally, O' is the grand-canonical potential of the refer- 
ence system. 

As a reference system for the VCA calculation, we 
take a cluster of four sites located in the (110) plane, 
and containing two from each inequivalent Cri and Cr2 
sites. [3^ We have checked that our results do not change 
substantially when using a cluster of 6 sites. We stress 
that the remaining hoppings, namely the ones between 
Cri-Cri along x and y are not neglected, but taken into 
account in the noninteracting lattice Green's function en- 
tering the Dyson equation in the expression for the VCA 
Green's function. Quite generally, multi-orbital strongly- 
correlated systems show a competition between different 



magnetic and orbital-ordered phases. In order to check 
that the ferromagnetic phase is the one with the low- 
est energy, we have compared its grand-canonical poten- 
tial with the one of different magnetic phases containing 
mixed antiferromagnetic and ferromagnetic components 
in different directions [IBl . Notice that the particle den- 
sity in the physical system is in general different from the 
one of the reference system ^] ; the former turns out to 
be slightly doped with n = 1.81 particles per unit cell. 
Moreover, for the parameters we have used here, which 
are the appropriate ones for Cr02, we confirm that the 
ferromagnetic state is more favorable energetically with 
respect to all antiferromagnetic states we have consid- 
ered, as well as to the paramagnetic state. 

In order to analyze the role of the different d or- 
bitals, we present in Fig. [T]the orbital and spin-resolved 
density of states (DOS). For comparison, we show in 
the inset the results obtained within the mean-field-like 
LDA-f U decoupling of the Hubbard interaction for the 
Cr-t2g orbitals. [Tsj In the LSDA-I-U, a pseudogap fea- 




FIG. 1: (Color online) Cr-t29 orbital-resolved density of states 
in the ferromagnetic state, calculated within the LDA-I-VCA 
approach for U=3eV, and J=0.9eV. NQP states are visible 
in the minority spin channel, just above the Fermi energy. 
The inset shows the LSDA-I-U density of states for the same 
parameters. 

ture opens at the Fermi energy in the majority (spin 
up) DOS, which is attributed to the crystal-field split- 
ting between the dyz+zx and dy^-zx orbitals. This pseu- 
dogap is already present in the (uncorrelated) LSDA re- 
sults shown in Fig. [21 although the effects of the inter- 
action in LSDA+U enhance this effect. On contrary, in 
the VCA calculation, correlation effects again reduce the 
pseudogap feature with respect to LSDA-I-U and shift it 
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to higher energies. 

A remarkable result of the VGA calculation is the pres- 
ence of a significant amount of density of states at the 
Fermi energy originating from the putatively localized 
dxy orbital. More specifically, we obtain that this or- 
bital is not completely half filled, its occupation being 
n^y « 0.87, while for the other two orbitals we have 
^y^+zx ^ 0.49, and n?'^"^^ !v 0.45. This is confirmed 
by our LSDA-I-DMFT calculation, discussed below, for 
which k:: 0.87, and n^^^^^ w 0.46. Consequently, the 
dxy orbital consists of quite itinerant electrons, although 
with a large effective mass, rather than of localized mo- 
ments. Notice that our findings are in contrast to pre- 
vious DMFT 17, [3] calculations, in which the Fermi 
energy only touches the high-energy tail of the dxy DOS, 
which can then be considered as localized moments. This 
is probably due to the large value of the interaction pa- 
rameter U used in these calculations. We believe that the 
smaller value of C/ ~ 3eV used here is more appropriate, 
as it is obtained from first principles. 

In the minority spin channel, NQP states are clearly 
visible predominantly in the dyz±zx orbitals just above 
the Fermi energy. The physical origin of NQP states is 
connected with the "spin-polaron" processes 2^, 21 1 : the 
spin-down low-energy electron excitations, which are for- 
bidden for half-metallic ferromagnets in the one-particle 
picture, turn out to be possible as superpositions of spin- 
up electron excitations and virtual magnons 23, 21 1 ■ An 
uniform superposition forms a state with the same total 
spin quantum number S = {N + l)/2 {N is the particle 
number in the ground state) as the low-energy spin-up 
state, but with one "spin-flip" , i. e. with z-component 
Sz — S — 1. If the Hamiltonian is spin- rotation invariant 
this state with one additional spin-down particle must 
have the same energy as the low-energy state with one 
additional spin-up particle, although its weight is reduced 
by a factor 1/N. We stress that spin-rotation invari- 
ance is crucial in order to obtain low-energy NQP states. 
For this reason, methods neglecting spin flip processes in 
Eq. ID) are not expected to provide a correct description 
of NQP states. Similar NQP states are obtained in the 
fully self-consistent LSDA-I-DMFT calculation described 
in the following paragraph. 

The above VGA calculation was carried out by start- 
ing from a model ([T]), whose parameters are obtained 
ab-initio from electronic-structure calculations. This pro- 
cedure does not take into account feedback effects of the 
correlation-induced charge redistribution into the LSDA 
potential. This is still a difficult task in a VGA calcu- 
lation. In addition, correlations and states in the Cg d 
orbitals are neglected. In order to explore the conse- 
quences of these effects, and to access finite temperatures, 
which is more difficult within a Lanczos-VGA approach, 
we also carried out a multi-orbital LSDA+DMFT cal- 
culation, which includes the full Gr d-manifold, as well 
as a complete spd-basis set for all the atoms in the unit 



cell. In addition, we include both self-energy and charge 
self-consistency, which means that many-body effects are 
taken into account in the evaluation of the LSDA poten- 
tials. Detailed of this method have been given elsewhere 
13 All previous LSDA-hU [H HI] or DMFT [B, [H 
studies, independently yield a narrow almost flat band of 
dxy character which produces the ferromagnetic phase of 
Gr02. We show that in contrast to these previous results, 
the fully self-consistent LSDA-I-DMFT results, in agree- 
ment with the above non-local VGA approach, yields an 
itinerant dxy orbital, which crosses the Fermi surface. It 
is interesting to note that despite the non-localized na- 
ture of the orbital, we still obtain a ferromagnetic phase. 
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FIG. 2: (Color online) Density of states obtained within the 
LSDA and LDSA-I-DMFT calculations for different tempera- 
tures. The inset shows the results for a smaller energy window 
around the Fermi level. 

The DOS obtained from our LSDA+DMFT calcula- 
tion is presented in Fig. [5] for two different values of 
T, and compared with the LSDA results. As discussed 
above, the LSDA Fermi level intersects the majority-spin 
bands near a local minimum and lies in the band gap of 
the minority spin states. Finite temperatures and corre- 
lation effects close this minimum around the Fermi level, 
as can be seen from the LDA-I-DMFT results in Figl^l 
No differences can be observed between the two DMFT 
results at different temperatures, except for the smear- 
ing of DOS features at higher temperature. For both spin 
channels, the DOS is shifted uniformly to lower energies 
in the energy range of -2 and -6eV, where predominantly 
the 0{p) bands are situated. This is due to the fact that 
correlated Cr{d) bands affect the 0{p) states through the 
Cr{d) — 0(p) hybridisation, so that the latter contribute 
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actively to the mechanism of the ferromagnetic ground 
state. The LSDA+DMFT calculation confirms the exis- 
tence of minority spin states just above the Fermi energy, 
as observed in the VGA calculation. 

We now show that many-body effects discussed above, 
especially the formation of NQP states, contribute signif- 
icantly to the energy dependence of the spin polarization 

P{E) = (N^iE) - N^iE))/iN^{E) + N^{E)) , (3) 

where N{E) is the density of states of majority f or mi- 
nority I electrons. Fig. [3] shows a comparison between 
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FIG. 3: (color online)Energy dependence of the spin polariza- 
tion obtained experimentally (Ref. [s^ decorated solid line) 
and by different theoretical calculations LDA (dot-dashed), 
DMFT (dashed), VGA (soUd) (Eq. ^. A broadening of 0.4 
eV, corresponding to the experimental resolution, has been 
added to the theoretical curves. 

the measured [39| and computed polarization for the dif- 
ferent ab-initio many-body calculations discussed in the 
present paper. Due to the tails of NQP states, polariza- 
tion is less than 100% even at the Fermi level. The LDA 
calculation clearly overestimates the spin polarization as 
it neglects correlation effects. On the other hand, both 
VGA and DMFT results show excellent agreement with 
the experiment at the Fermi level. At an energy of about 
0.5 eV polarization is reduced by about 50%, and up to 
this energy the agreement of the VGA calculation is ex- 
cellent, while DMFT results overestimate depolarisation 
effects away from the Fermi energy. The disagreement 
for energies above « 0.5 eV is probably due to terms not 
included in the VGA Hamiltonian ^ , such as Cg orbitals 
which start becoming important at higher energies. 

In conclusion, we addressed the effects of electronic 
correlations in Cr02 and showed how these consider- 
ably change the picture obtained within LSDA and 



LSDA-f U, despite the fact that the interaction is not too 
strong. More specifically, while in LSDA-I-U exchange 
and crystal-field splitting pin the Crd^y electrons to be- 
come localized moments, the competition with charge 
fluctuation effects, which are taken into account in our 
DMFT and VGA calculations, induce a ferromagnetic 
state in which also dxy electrons are itinerant although 
with a large effective mass. We believe that this is the 
correct behavior at this intermediate value of the interac- 
tion parameters, which we derive ab initio, in contrast to 
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results obtained with a larger value of U 
We show that correlation-induced NQP states as well as 
an appropriate description of their nonlocal contributions 
to the self energy as obtained by VGA are crucial for a 
correct description of the energy dependence of the po- 
larization. Additional effects, such as, e. g., disorder 
or phonons, are also expected to contribute to the spin 
depolarisation. 
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